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Abstract 

The precision and discovery potential of a neutrino factory based on muon 
storage rings is studied. For three-family neutrino oscillations, we analyse how 
to measure or severely constraint the angle ^13, CP violation, MSW effects and 
the sign of the atmospheric mass difference Arn^^. We present a simple ana- 
lytical formula for the oscillation probabilities in matter, with all neutrino mass 
differences non-vanishing, which clarifies the subtleties involved in disentangling 
the unknown parameters. The appearance of "wrong-sign muons" at three ref- 
erence baselines is considered: 732 km, 3500 km, and 7332 km. We exploit the 
dependence of the signal on the neutrino energy, and include as well realistic 
background estimations and detection efficiencies. The optimal baseline turns 
out to be 0(3000 km). Analyses combining the information from different base- 
lines are also presented. 
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1 Introduction 



The atmospheric |jT|, plus solar |^ neutrino data point to neutrino oscillations 0, |^ 
and can be easily accommodated in a three-family mixing scenario. 



Let U, with (i/g, t'u, '^r)"^ = U ■ (z/i, z/2, be the leptonic Cabibbo-Kobayashi- 



Maskawa (CKM) matrix in its most conventional parametrization 




U = U23Ui3Ui2 = c, s,, I I 1 



/ C12 S12 o~ 
-S12 C12 



ci3 y V 1 



with S12 = sin ^12, and similarly for the other sines and cosines. Oscillation experiments 
are sensitive to the neutrino mass differences and the four parameters in the mixing 
matrix of Eq. (|l|): three angles and the Dirac CP-odd phase. 

The SuperKamiokande data on atmospheric neutrinos are interpreted as oscil- 
lations of muon neutrinos into neutrinos that are not z/g's, with a mass gap that we 
denoteQby Am^s. Roughly speaking, the measured mixing angle ^23 is close to maxi- 
mal and lAmlgl is in the range 10^^-10"^ eV^. The solar neutrino deficit is interpreted 
either as MSW (matter enhanced) oscillations or as vacuum oscillations (VO) 0] 
that deplete the original u^s, presumably in favour of z/^'s or alternatively into ster- 
ile neutrinos. The corresponding squared mass differences -(9(10~^-10~^) eV^ for the 
large mixing angle MSW solution (LMA-MSW), O(10~^) eV^ for the small mixing 
angle MSW solution (SMA-MSW), or C(10"^°) eV^ for VO- are significantly below 
the range deduced from atmospheric observations. We identify this mass difference 
with Am^2 this parametrization. Its sign is constrained by solar data: while the 
SMA-MSW solution exists only for positive Am^2) the LMA-MSW range there is 
also a small window at negative values 

These oscillation signals will be confirmed and further constrained in ongoing and 
planned atmospheric, solar and long baseline reactor experiments |^, as well as in 
future long baseline accelerator neutrino experiments 0. In a few years they will 
answer the question of sterile neutrinos contributing or not to present data. The MSW 
effect is expected to play a major role in explaining the solar deficit and both solar and 
reactor experiments will also clarify whether Nature has chosen the LMA-MSW rather 
than SMA-MSW or VO solutions. 

The atmospheric neutrino parameters will be known with better precision as well. 
Experimental information relevant for a more precise knowledge of the atmospheric 
neutrino fluxes will be available |l^, [ill. Also, projected long baseline accelerator 



experiments will improve the precision of lAmggI and 6*23. For instance, |Am23| is 
expected to be measured at MINOS with an accuracy below 10% if lAm^gl > 3 x 10~^ 
eV2 0. 



^ Amfj = nij — raf throughout the paper. 
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Nevertheless, there is a strong case for going further in the fundamental quest of the 
neutrino masses and mixing angles, as a necessary step to unravel the fundamental new 
scale(s) behind neutrino oscillations. In ten years from now no significant improvement 
is expected in the knowledge of: 



The angle 613, which is the key between the atmospheric and solar neutrino 
realms, for which the present CHOOZ bound is sin^ ^13 < 5 x 10~^ |T^[Q. 



The sign of Am"^^, which determines whether the three-family neutrino spectrum 
is of the "hierarchical" or "degenerate" type (i.e. only one heavy state and two 
almost degenerate light ones, or the reverse). 

Leptonic CP-violation. 

The precise study of matter effects in the u propagation through the Earth: 
a model-independent experimental confirmation of the MSW effect will not be 
available. 



The most sensitive method to study these topics is to measure the transition prob- 
abilities involving z/g and z/g, in particular i^e('^e) This is precisely the golden 
measurement at the neutrino factory [|l^. Such a facility is unique in providing high 
energy and intense t'e(^'e) beams coming from positive (negative) muons which decay 
in the straight sections of a muon storage ring |]16[. Since these beams contain also 
i^fii^fj.) (but no i/^(z/^)!), the transitions of interest can be measured by searching for 
"wrong-sign" muons: negative (positive) muons appearing in a massive detector with 
good muon charge identification capabilities. 

The first exploratory studies of the use of a neutrino beam with these characteristics 



T7| were done in the context of two-family mixing. In this approximation, the wrong- 



sign muon signal in the atmospheric range is absent, since the atmospheric oscillation 
is ^ Vt-. The enormous physics reach of such signals in the context of three-family 

where the authors put the emphasis on the 



neutrino mixing was first realized in 
measurement of the angle ^13 and CP- violation (see also |19[)- The latter may be at 
reach if the solar deficit corresponds to the LMA-MSW solution |ll8| , ^ |2l|]. Recently, 
it has also been shown p3 that the precision in the knowledge of the atmospheric 



parameters 6*23 and 



|Am^3| 



can reach the percent level at a neutrino factory, using 



muon disappearance measurements. Furthermore, it was pointed out p2|. 



that 



the sign of Am23 be determined at long baselines, through sizeable matter 

effects. The importance of measuring precisely Earth matter effects in a clean and 
model-independent way, both for the understanding of the fundamental parameters 
and for their intrinsic interest has been stressed in refs. 12^, |2B|, JTt 



The aim of this paper is to identify the optimal baselines for studying the above 
itemized topics. This requires to include in the analysis the maximum information that 
can be attained at any fixed baseline. While all previous analyses have been based in 
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energy-integrated quantities, we will take into account the neutrino energy dependence 
of the wrong-sign muon signals, together with the information obtained from running 
in the two different beam polarities. 

We will consider in turn scenarios in which the solar oscillation lies in the SMA- 
MSW or VO range and in the LMA-MSW range. In the latter, the dependence of the 
oscillation probabilities on the solar parameters 6'i2, Am^2 5 the CP-odd phase, 

S, is sizeable at terrestrial distances and complicates the measurement of 613 due to the 
presence of other unknowns (mainly 6). This potential difficulty was first pointed out 
in 



20| . The previous analysis of the sensitivity to ^13 neglected solar parameters 



and is thus only valid for the SMA-MSW and VO solutions, or the LMA-MSW if 613 is 
large enough. In the present paper a higher statistics is considered, allowing to explore 
smaller values of di^, and the remark is very pertinent. We will discuss in detail the 
issue of how to disentangle 6*13 and 6, guided by an approximate analytical formula for 
the oscillation probabilities in matter including two distinct mass differences. As we 
will see, the choice of the correct baseline is essential to solve this problem. 

We shall consider the following "reference set-up": neutrino beams resulting from 
the decay of 2 x 10^°/i"'"'s and/or /i~'s per year in a straight section of an = 50 
GeV muon accumulator. A long baseline (LBL) experiment with a 40 kT detector 
and five years of data taking for each polarity is considered. Alternatively, the same 
results could be obtained in one year of running for the higher intensity option of 
the machine, providing 10^^ useful /i"'"'s and /i~'s per year. A realistic detector of 
magnetized iron will be considered and detailed estimates of the corresponding 
expected backgrounds and efficiencies included in the analysis. Three reference detector 
distances are discussed: 732 km, 3500 km and 7332 km. 



A preliminary version of this work was presented in |30] 



2 Fluxes and charged currents. 
2.1 The Neutrino Factory 

One of the most encouraging outcomes of the recent neutrino factory workshop at 



Lyon RT| was the convergence of the various machine designs existing previously, to an 



essentially unified design |^, based on a muon accumulator with either a triangular 



or a bow-tie shape. Both geometries permit two straight sections pointing in different 
directions, allowing two different baselines. 

It was also agreed in Lyon that the beam power on target should not exceed 4 MW. 
This, in turn, limits the production of muons to 10^^ per year, out of which only 20- 
25 % are useful (that is, decay in the straight sections pointing towards the detectors). 
Agreement was also found upon other important parameters, namely the machine dual 
polarity (i.e, the ability to store and although not simultaneously) and the 
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fo{x) 






2x^(3 -2x) 


2x^(1 - 2x) 




12x2(1 - x) 


12x^(1 - x) 



Table 1: Flux functions in the muon rest- frame as in ref. f^3^. 

maximum realistic energy to accelerate the stored muons, which was fixed at 50 GeV. 

Ultimate sensitivity to the neutrino mixing matrix parameters, in particular to 
and (5, require a data set as large as possible. The analysis presented in this paper 
assumes a total data set of 10^^ useful /x"^ decays and 10^^ useful /i~ decays. 



2.2 Number of events 



In the muon rest-frame, the distribution of muon antineutrinos (neutrinos) and electron 
neutrinos (antineutrinos) in the decay /i^ — >• + feipe) + ^niy^) is given by: 



dxdfl 



— [fo{x)TVMx)cose] 



(2) 



where denotes the neutrino energy, x = 2Ei^/m^ and is the average muon 
polarization along the beam directions. 6 is the angle between the neutrino momentum 
vector and the muon spin direction and is the muon mass. The positron (electron) 
fiux is identical to that for muon neutrinos (antineutrinos), when the electron mass is 
neglected. The functions /o and /i are given in Table |I] PB| . 

In the laboratory frame, the neutrino fiuxes, boosted along the muon momentum 
vector, are given by: 



dydQ 



dydfl 



4:n,. 



Ey (1 



(3 cos if) I 



P cos (f) 



TV^[ml -AEly {1-/3 COS ip)]} , 
2^^" ^ ^ l-2Ely{l-Pcosip) 



Ey (1 



TV, 



m, 



(3 cos if) I 
2Ely {I -13 cos ^)\} . 



(3) 



Here, (3 = -^1 — m?^/ EJ^, E, is the parent muon energy, y = E^/E,, is the number 
of useful muons per year obtained from the storage ring and L is the distance to the 
detector, ip is the angle between the beam axis and the direction pointing towards 
the detector, assumed to be located in the forward direction of the muon beam. The 
angular divergence will be taken as constant, S^p ~ 0.1 mr. 

Unlike traditional neutrino beams obtained from tt and K decays, the fiuxes in 
Eq. (^, in the forward direction, present a leading quadratic dependence on Ey. As a 
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consequence, the oscillation signal does not decrease with increasing E^. In Appendix 
A we present our numerical results for the i^e(^'e) and i^^ii^fj,) fluxes. 

The charged current neutrino and antineutrino interaction rates can be computed 
using the approximate expressions for the neutrino-nucleon cross sections with an 
isoscalar target, 

E E 
a^N ~ 0.67 X 10"^^ X ^ X m^ , a^^N ~ 0.34 x lO'^^ ^ _^ ^ ^2 

It follows that the number of charged current (CC) events at a neutrino factory scales 
cubically with energy. In Appendix A we also include our numerical results for the 
rates of e"^ and fi"^ production, from a fi^ beam. 



3 Oscillation Probabilities 
3.1 In vacuum 

Atmospheric or terrestrial experiments have an energy range such that Am^ L/ Ey 
for the smaller (Am^j) but not necessarily for the larger (Amjs) of these mass gaps. 
Even then, solar and atmospheric transitions are not (provided ^13 7^ 0) two separate 
two-generation oscillations. For lAm^gl |Am23|, neutrino oscillation probabilities at 
terrestrial distances are accurately described by only three parameters, ^23? Am23 = 
Am^3 and 6'i3: 



Py.u^iueu^) = ^23 sin 2^13 sin' 
Pu,ur(u,Pr) = C23 sin^ 26*13 sin^ 



Ai3 


L 






Ai3 


L 


2 




Ai3 


L 


2 



where 



'u^Mu^u^) = c^3 sin^ 2^23 sin^ ( ) , (5) 

A- = ^ (6) 

Eqs. d^) are a very good approximation when the solar parameters lie in the SMA- 
MSW or VO range. The present best flt value for ^13 is in the range 6°-8°p^, |35 



although it is compatible with zero within errors. Among the transitions in Eq. (J^) the 
channels Ue v^iI't have clearly the best sensitivity to a small ^13. Experimentally, 
the measurement of z/g — > z/^ oscillations through the appearance of wrong-sign muons 



is far superior to that of Ve oscillations through r detection. In it was 

shown that the sensitivity to ^13 of the former can improve the present limits, which 
are mainly set by Chooz |13|, by at least two orders of magnitude. 

At the neutrino factory, precision measurements for |Am23| and 6*23 can also be 
performed. Measurements of r appearance or /j, disappearance may be competitive 
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with wrong-sign /i signals for small values of ^^13, because of the cosine dependence of 
the corresponding probabilities in Eqs. We do not develop this topic further in the 
present work, see . 

In the LMA-MSW scenario, for which the effects of Am^2 may be relevant at the 
neutrino factory, a good and simple approximation for the Uf, transition proba- 

bility is obtained by expanding to second order in the small parameters, 6'i3, A12/A13 
and A12L: 

22 2 f A13 L\ 2 2 if A12 L 

Pu,u^{Peu^) = ^23 sin 2^13 sin ^ ) + C23 sin 2^12 sin y-^— 

+ J cos (±5 - sm [-^) , (7) 

where here and throughout the paper the upper/lower sign in the formulae refers to 
neutrinos/antineutrinos, and 

J = Ci3 sin 26*12 sin 26*23 sin 26*13 (8) 

is the usual combination of mixing angles appearing in the Jarlskog determinant. The 
first term in Eq. (|^ is quadratic in sin 6*13, whereas the leading term in A12 is linear. 



The latter may then be significant or even dominant for very small values of ^13 |20 
when Ami2 and 6*12 are in the range allowed by the LMA-MSW solution. 

At "short" distances, such as 732 km, Eq. (|^ can be further approximated by: 
Pu^Mu^u,) = S23 sm 2^13 — — +J cos5 — — +C23sm 2^12 — — . (9) 



The CP-odd term in the probability (i.e. the one proportional to sin 5) has dropped 
because it is of higher order in Aij L. The comparison of the two polarities is then not 
useful (except for doubling the statistics ), since the CP-conjugated channels measure 
the same probability. Furthermore, all terms in Eq. (^ have the same dependence 
on the neutrino energy and the baseline, and consequently is very hard to disentangle 
them. As a result we expect a large correlation between the parameters 6*13, 6, A12. 
Even though long baseline (LBL) reactor experiments will provide a measurement 
of |Am^2l if the LMA-MSW solution is at work, the parameters 6*13 and 6 will have to 
be determined simultaneously at the neutrino factory. It will then be necessary to go to 
longer baselines where the energy dependence of the different terms in Eq. (^ differs, 
and where the comparison of the neutrino and antineutrino probabilities provides non- 
trivial information to separate 6 from ^13. 

It is uncertain 0] whether solar experiments will determine the sign of Am^2 if it 
lies in the LMA-MSW range. If it remains unknown, it implies an ambiguity in the 
determination of 6: to reverse the sign of Ami2 in Eq. (0) is equivalent to replacing 6 
by 5 -|- TT. Notice, though, that whether there is CP-violation or not is independent of 
whether the phase in any parametrization is 5 or 5 + tt. 
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It is possible to construct a measurable CP-odd asymmetry, which in vacuum is pro- 
portional to sin 5. In refs. [|1^ and ||21|, the authors considered the following integrated 
asymmetry (see also p7| ): 




{N[f,-]/N,[e-]U-{N[f^+]/N^[e+]} 
{iV[/i-]/iV,[e-]}+ + {iV[/i+]/iV,[e+]} 



(10) 



The sign of the decaying muons is indicated by a subindex, N[^^]{N[^^]) are the 
measured number of wrong-sign muons, and iVo[e+](Ai'o[e~]) are the expected number 
of Uei^e) charged current interactions in the absence of oscillations. The significance 
of this asymmetry (i.e. the asymmetry divided by its statistical error) scales with the 
baseline and neutrino energy in the following way: 



The best sensitivity to a non-zero CP-odd asymmetry is found at the maximum of 
the atmospheric oscillation. At the corresponding distance, however, matter effects are 
already important and should be taken into account, as we proceed to discuss in the 
next subsection. 

3.2 In matter 

Of all neutrino species, only Ue and z/g have charged-current elastic scattering ampli- 
tudes on electrons. This, as is well known, induces effective "masses" fj, = ±2 Ei,A, 
where the signs refer to and and A is the matter parameter. 



with rie the ambient electron number density P]. 

Matter effects may be important if A is comparable to, or bigger than, the quantity 
Aij for some mass difference and neutrino energy, and if distances are large enough for 
the probabilities to be in the non-linear region of the oscillation. 

For the Earth's crust, with density p ~ 2.8 g/cm^ and roughly equal numbers of 
protons, neutrons and electrons, A ~ 10"^'^ eV. The typical neutrino energies we are 
considering are tens of GeVs. For instance, for E^, = 30 GeV (the average Ue energy 
in the decay of E^^ = 50 GeV muons) A = 1.1 x 10^^ eV^/GeV ~ A23. This means 
that matter effects will be important at long distances. Notice that at L = 732 km 
and 3500 km the neutrino path remains in the Earth crust, whereas for 7332 km the 
deeper flight path meets a denser medium and A = 1.5 x 10~^ eV^/GeV. 




(11) 



A = V2GFne, 



(12) 
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Figure 1: The signal-over-noise ratio of the CP-odd asymmetry of Eq. ^T^ ) as a function 
of the distance. The thick line corresponds to the integrated asymmetry, while the 
dashed lines correspond to the asymmetry computed in five energy bins of equal width 
AE,, = 10 GeV. The neutrino mixing parameters correspond to the S MA- MS W solution 
to the solar anomaly: Am^g = 2.8 x 10"^ eV^, Am^2 = 6 x 10~^ eV^ , sin^ 26^12 = 
6 X 10-^ = 45°, = 13° and 5 = 90°. The muon energy zs E^ = 50 GeV and the 
matter parameter A is varied with the distance as in [S^j. 



3.2.1 Neglecting solar parameters: VO or SMA-MSW solutions 

Consider the case when A 12 is negligible compared to A23, and A. In the approx- 
imation of constant rie, the transition probability in matter governing the appearance 
of wrong-sign muons can then be read from p8|] : 

^....(PeP.) ^ ^23 sin^ 2^13 sin^ (^^^ , (13) 

where 

= v/[Ai3 cos 2^13 TAf+ [Ai3 sin 2^^i3]' • (14) 

For A = 0, Eq. (|T^) reduces to the corresponding vacuum result: the first line in 
Eq. (H). At short distances, that is for BzfL/2 sufficiently small, the sinus in Eq. ( jlB]) 
can be expanded and 

/A Zy\ ^ 

^^e^M(PePM) ~ ^23 SIH^ 2^13 j > (15) 

which also coincides with the vacuum behaviour for small A13L/2, even when A ^ 
Ai3Cos26'i3. Eq. (plSj) is a good approximation up to ~ 3000 km [1^. Note that the 
leading dependence on 6^13 is quadratic as in vacuum. 
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In contrast with the vacuum result, the probabihty in matter depends on the sign of 
Amfg. It follows from Eqs. (0) and (|14D that a change on this sign is equivalent to a CP 
transformation^, that is, interchanging the probability of neutrinos and antineutrinos. 
Thus matter effects induce by themselves a non-vanishing CP-odd asymmetry, and 
the best sensitivity to the sign is achieved when the sensititivity to the asymmetry of 
Eq. (p!OD is maximal. Fig. || shows the significance of this asymmetry as a function of 
the baseline. The maximum sensitivity to the sign is thus expected at (9(7000) km, 
although it is already very good at much shorter distances: notice the large number of 
standard deviations in the y axis. 

It is important to stress that having sizeable matter effects does not necessarily 
imply having sensitivity to the sign. For example if A 3> Ai3Cos26'i3, the sensitivity 
is lost, even though matter effects are important at large distances. The optimal 
sensitivity occurs for A ~ A13 cos 2^13, which is an energy dependent condition. In 
Fig. ID the asymmetry resulting from each of five energy bins of width AEy = 10 
GeV is also shown: the asymmetries in different bins peak at different distances. This 
dependence suggests that using the information in energy bins can further improve the 
measurement of the sign of Am23. 



3.2.2 With solar parameters: LMA-MSW solution 

In the LMA-MSW solar scenario, the effects of Amfg not negligible over terrestrial 
distances, given the high intensity of the neutrino factory. The exact oscillation proba- 
bilities in matter when no mass difference is neglected have been derived analytically in 
40| . However, the physical implications of the formulae in |^ are not easily derived. 



A convenient and precise approximation is obtained by expanding to second order in 
the following small parameters: 6^13, A12/A23, A12/A and A12 L. The result is (details 
of the calculation can be found in appendix C): 

Pueu^{v,u^) = S23sm 2^13 l-g-l sm |— ^1 + C23 sm 2^12 j sm 

^Ai2Ai3 . (AL\ . (B^L\ / A13L 



^ '^t^''^ [-) - [-r ) i^* - • 

where = \A ^ ^isl- Once again, this expression reduces to the vacuum result, 
Eq. d^), in the limit A ^ 0. As already remarked in the previous subsection, a reversal 
of the sign of A 12 in Eq. (piBf ) can be simply traded by an shift of tt in 5, and we will 
stick to positive A12 in the numerical exercises below. 



We have numerically compared Eq. (0) with the exact formulae of ^y], in the 
range 1° < ^13 < 10°. Consider for instance the average energy = 30 GeV and the 
following set of values: Amj^ = 2.8 x lO'^ eV^, Ou = 22.5°, 623 = 45° and 6 = 90°. For 



^Recall that in this approximation A 
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Figure 2: a) Scaling with L of the ratio of the number of wrong-sign muons fi" induced 
solely by P2 in Eq. ( fT^ j to the statistical error in the total number of wrong- sign muons, 
for 5 = 0°; b) the same for wrong-sign muons fi~ induced by P3, for 6 = 90°. The 
normalization of the y axis is arbitrary. The oscillation parameters are: Am"^^ = 
2.8 X 10-3 eV^, Amj^ = 1 x 10"^ eV^, Ou = 22.5°, 623 = 45° and ^13 = 8°. 



Aml^ = 1 X 10"^ eV2, the difference is < 10 % (< 20 %) at L ~ 3000 (7000) km. For 
Amj^ = 1 X 10"^ eV^, the error diminishes to < 2.5 % (< 10 %) at L ~ 3000 (7000) 
km. Shghtly better accuracy is obtained for 6 = 0°. 

As before, matter effects in Eq. (|T^) induce an asymmetry between neutrinos and 
antineutrinos oscillation probabilities even for vanishing 6. For this reason the CP-odd 
asymmetry, Eq. ([To|) , is not the most transparent observable to determine the optimal 
distance for measuring 6. A better way of addressing the issue is to ask at what distance 
the significance of the terms which depend on S is maximal. 

For 6*13 > 1°, the dominant contributions in Eq. (|16|) are: 
Pi = 4 sin^ 2^13(^1 

P2 ^ ^— ^-os5cos(^^j ^^H— r^i^j' 

_ ^Ai2Ai3 . . . (A^3L\ . (AL\ . (B^L\ 
P3 = ±J — -g- sm5sm (^^j sm j sm j . (17) 

In Fig. 1^ we show the significance of the 5-dependent terms, P2 and P3, as function 
of L. The significance is defined as the fraction of wrong-sign muons-for positive muons 
in the beam- resulting from a given term, over the statistical error in the measurement 




10 




0^ 

X 



0.4 



0.3 



0.2 



0.1 



V 



10 



20 



(a) 



30 
, [GeV] 

(b) 



40 



50 



Figure 3: Neutrino energy dependence of the different terms in Eqs. p^j.- E^Pi (solid 
line), EIP2 6 = 0° (dashed Ime) for and El P3 for 5 = 90° (dotted Ime) , at L = 132 
km, for neutrinos (a) and antineutrinos (b). The parameters correspond to the LMA- 



MSW solution: Am^^ 



2.8 X 10-3 eV^, A 



m 



12 



1 X 10" 



12 



22.5°,^ 



23 



45= 



and 9 



13 



X 



ttf 2 




20 30 40 

[GeV] 

(a) 



[GeV] 
(b) 



Figure 4: The same as in Fig. ^ but for L = 3500 km. 



of the total number of wrong-sign muons. Results for several values of the neutrino 
energy are depicted. 

The term in cos 5, P2, is more significant than P3 at short distances. Unfortunately, 
this sensitivity to 5 through P2 is fake, because at short distances there is no way to 
separate P2 from the leading term. Pi: they have similar energy dependence and do 
not differ in the two polarities, as illustrated in Fig. |^. In order words, a change in 5 
can be compensated by a change in 6*13 to keep the total probability unchanged. 
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At 3500 km the situation is very different, though: as Fig. ^ illustrates, the energy 
dependence of the three terms is quite different and furthermore P3 -which changes 
sign with the beam polarity- is considerably larger. The comparison of the number of 
wrong-sign muons detected running in the two polarities and the binning in energy of 
the signal are thus strong analysis tools to disentangle ^13 and 6 at baselines around 
3500 km. Notice that this optimal distance is in nice agreement with the previous 
studies based on the significance of integrated asymmetries [jl8|, ^ updated in 
Fig. ^ for the present set-up. 

The summary of this long discussion is that baselines much larger than 732 km are 
needed, for the following reasons: 



• In the SMA-MSW or VO scenarios, the sign of Aml^ can only be determined 
for distances such that matter effects are sizeable and the CP asymmetries they 
induce measurable. This happens aX L = O (3000 km) or larger. 

• In the LMA-MSW scenario, there is a strong correlation between ^13 and 6 at 
short distances. It is necessary to go far so that the terms most sensitive to them 
show a different energy dependence, and the signals in the CP-conjugate channels 
differ sizeably, allowing the simultaneous measurement of both parameters. 

These expectations will be sustained by a detailed energy analysis in the following 
sections. 



4 Detection of wrong sign muons 
4.1 A Large Magnetized Calorimeter 

For the present study we consider a Large Magnetized Iron Calorimeter such as the 
one proposed in [Q. The apparatus, shown in Fig. |^ is a huge cylinder, of 10 m 
radius and 20 m length. It is made of 6 cm thick iron rods intersped with 2 cm thick 
scintillators segmented along their length. Its fiducial mass is 40 kT. A superconducting 
coil generates a solenoidal magnetic field of 1 T. 

The detector axis is oriented to form a few degrees with the direction of the neutrino 
beam. Thus, a neutrino crossing the detector sees a sandwich of iron and scintillator. 
The coordinates of the track transverse to the cylinder axis are measured from the 
location of the scintillator rods, while the longitudinal coordinate is measured from 
their segmentation. As discussed in [^, the performance of the device would be 



similar to the one expected for the MINOS detector P2 |. 

Neutrino interactions in such a detector have a clear signature. For example, a z/^ 
charged current (CC) event is characterized by a muon, which is seen as a penetrating 
long track, plus a shower resulting from the interactions of the hadrons in the event. 
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Figure 5: The signal- over-noise ratio of the CP-odd asymmetry of Eq. as a func- 
tion of the distance, after subtracting the fake matter induced asymmetry. The thick 
line corresponds to the integrated asymmetry, while the dashed lines correspond to the 
asymmetry computed in five energy bins of equal width AEi, = 10 GeV. The neu- 
trino mixing parameters correspond to the LMA-MSW solution to the solar anomaly: 
Amis = 2.8 X 10-3 eV^, Amj^ = 1 x 10"^ eV^, 6^ = 22.5°, ^23 = 45°, ^13 = 13° and 
6 = 90°. The muon energy is = 50 GeV and the matter parameter A is varied with 
the distance as in ^3^ . 




Figure 6: Sketch of the Large Calorimeter for the Neutrino Factory. 

Fitting the muon track determines its charge and momentum, and a calorimetric mea- 
surement allows the determination of the hadronic momentum vector. In contrast, 
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for a z/g CC event the electromagnetic shower due to the prompt electron cannot be 
disentangled from the hadronic shower in an event-by-event basis, and therefore, the 
event looks similar to a neutral current (NC), which is characterized by having no 
penetrating track. 

A realistic simulation of the response of an iron calorimeter must be able to compute: 
1) whether the primary muon characterizing a CC event was identified or not, 2) 
whether non prompt muons arising from the decays of hadrons were missidentified as 
primary muons and 3) the muon and hadronic momentum vectors. 

To address the above points we have written a Monte Carlo simulation based in the 
GEANT 3 package [H3|. The apparatus is simulated with the correct geometry, and 



neutrino interactions are generated at random points in the fiducial volume. Then, 
every particle produced in the interaction is followed until it decays, exits the detector 
or undergoes a nuclear interaction. 

When the muon track is very short, it cannot be disentangled from the other tracks 
in the hadronic jets. The peak of the hadron shower occurs at about 10 cm from the 
interaction vertex and essentially no hadronic activity remains at 100 cm from it. We 
thus impose the conservative criteria that, in order to be reconstructed, a muon track 
must be longer than 100 cm. About 99.2 % of all z/^ CC events at 50 GeV produce 
primary muons that satisfy this condition. 

All muons, either primary or arising from the decay of hadrons, are tracked through 
the entire volume, and a hit is recorded each time a scintillator rod is crossed. The 
tracks can then be fitted to obtain the muon momentum. However, once the distance 
travelled by the muon is known, it is more efficient to use a simple smearing, which 
takes into account correctly both the effect of detector resolution and the multiple 
scattering. 

As an illustration of the smearing procedure, the relative error in the measurement 
of the muon momentum as a function of the muon momentum is shown in Fig. |^. 
The resolution decreases rapidly for low momentum muons (for < 5 GeV, a much 
better resolution would be obtained from the measurement of the muon range), and 
improves rather smoothly for large momentum {5P^/P^ ~ 4% for ~ 7 GeV, while 
5P^/P^ ~ 3% for P^ ~ 50 GeV). The fact that (5P^/P^ is almost constant for large 
momentum is due to the dominance of the multiple scattering term in the resolution. 

Hadrons, in contrast, are followed until they decay or undergo a nuclear interaction. 
In the latter case their momenta are added to the hadronic energy vector. Finally, 
both the magnitude and direction of the hadronic vector was smeared to account for 
the resolution of the detector. For details we refer to . 



A large sample of neutrino interactions, corresponding to 10 charged and neutral 
currents, 5 x 10^ z/g charged and neutral currents and the same data samples for the 
opposite polarity were analyzed for this study. 
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Figure 7: Relative error in the measurement of the muon momentum as a function of 
the muon momentum. 



4.2 Search for wrong-sign muon events 



We consider first the {u^ , z/g ) neutrino beams originating from a /x^ beam of = 50 
GeV (the dependence of the backgrounds with the energy of the muon beam was 
discussed in where it was shown that optimal performance is obtained at the 
highest possible energy). The bulk of the events in the detector are 77^ charged currents, 
signaled by the presence of a positive primary muon in the event, I^^ and z/g neutral 
currents, which are events with no primary leptons, and z/g charged currents, for which 
we assume that the primary electron is not identified. On top of those events, one 
searches for wrong sign /i~ arising from the z/^ produced via the oscillation z/g —>■ z/^. 
Table ^ shows the number of interactions corresponding to a total of 10^^ useful /i"*" 
decays and a 40 kT detector at our reference baselines. For illustration, we will consider 
the oscillation parameters for the signal in the LMA-MSW scenario: Am23 = 4 x 
lO-^eV^, Amj^ = 10"^ eV^, ^13 = 13°, ^12 = 22.5° and 623 = 45°, as in the tables 
of Appendix B. Notice however that our results for the fractional background and 
efficiencies should be rather insensitive to the particular choice of parameters. 



Baseline (km) 


CC 


z/g CC 


17^ + Z/g NC 


z/^ (signal) 


732 


3.5 X 10^ 


5.9 X 10^ 


3.1 X 10^ 


1.1 X 10^ 


3500 


1.5 X 10^ 


2.6 X 10^ 


1.3 X 10*^ 


1.0 X 10^ 


7332 


3.5 X 10^ 


5.9 X 10^ 


3.0 X 10^ 


3.8 X 10^ 



Table 2: Data samples expected in a 40 kT detector for 10^^ useful jj,'^ decays, z/g — z/^ 
oscillations with parameters as in Appendix B. 

The potential backgrounds to the wrong sign ^~ events signaling the presence of 
oscillations are: 
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Figure 8: Definition of the kinematical variables used in this study. 

1. 17^ CC events in which the positive muon is not detected, and a secondary negative 
muon arising from the decay of tt~,K~ and D~ hadrons fakes the signaL The 
most dangerous events are those with D~ — > ^~ , which yield an energetic muon 
with a spectrum similar to the signal. 

2. i/g CC events, for which it is assumed that the primary electron is never de- 
tected. Charm production is not relevant for this type of events since the charmed 
hadrons in the hadronic jet are predominantly positive. Instead, fake fi~ arise 
from the decay of negative pions and kaons in the hadronic jet. 

3. 77^ and z/g NC events. Fake /i~ arise in this case also predominantly from the 
decay of negative pions and kaons, since charm production is suppressed with 
respect to the case of CC. 

At first sight these backgrounds seem impressive. Fortunately, simple kinematical 
cuts can suppress them very efficiently. One exploits the fact that for signal events the 
/i~ candidate is harder and more isolated from the hadronic jet than for background 
events. We thus perform a simple analysis based in two variables, namely: the momen- 
tum of the muon, P^, and a variable measuring the isolation of the muon, Qt = Pfj, sin 9 
(see Fig. ||). 

To illustrate the rejection power of this analysis. Fig. |^ shows the efficiency for 
signal detection and the fractional backgrounds as a function of and Qt for z7^ 
charged and neutral currents. Also shown is the signal-to-noise ratio, S/N, defined 
as the ratio of the signal selection efficiency and the error in the subtraction of the 
number of background events that pass the cuts, Nt. The error is taken to be Gaussian 
for large Nf, (~ and Poisson otherwise. Notice that charm production is the 

dominant background for z7^j CC, while vr decay dominates the NC backgrounds. 

Inspection of Fig. |^ shows that the S/N is rather flat for > 5 GeV and Qt > 
0.5 GeV. Cutting at P^ > 7.5 GeV, Qt > 1.0 GeV, maximizes S/N. However, given 
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z/^Charged Currents 




Figure 9: Signal and backgrounds for z/^j CC and NC events. 
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the flatness of tlie S/N ratio, one can vary tliese values generously with very little 
difference in the final results. Table ^ shows the fractional background contamination 
and the signal efficiency, while Table § shows the number of background and signal 
events that pass the cuts for the data sample discussed above. Notice that the residual 
backgrounds are quite sizeable at L = 732 km, small at L = 3500 km and negligible at 
L = 7332 km. It is possible to optimize the analysis for very long baseline in order to 
achieve higher efficiency. However, we have chosen to use the same set of cuts for the 
three baselines. 





I/e CC 


T7^+ Z/eNC 


//^(signal) 


1.0 X 10-^ 


5.0 X 10-^ 


1.0 X 10-^ 


3 X 10-1 



Table 3: Fractional backgrounds and signal selection efficiency for the wrong sign muon 
search with decays. 



Baseline (km) 


CC 


CC 




i/^ (signal) 


732 


3.5 


30 


31 


3.3 X 10^ 


3500 


0.1 


1.2 


1.2 


3 X 10^ 


7332 


< 0.1 


< 0.2 


< 0.2 


1.3 X 10^ 



Table 4: Events surviving the cuts in a 40 kT detector for 10^^ useful decays. 
Oscillation parameters as in Appendix B. 



A potential source of fake wrong sign muons not discussed here is that due to a 
wrong measurement of the charge. In it was estimated that, for = 50 GeV, 
this background could be reduced to a very small level, of the order of 10"^ or less. We 
have not included this source of background in the analysis. 

The same exercise has to be repeated when a fi~ beam is considered. The resulting 
neutrino beams are now z/^ and z/g and the signal events are z7^ . Fig. ^ shows the 
efficiency for signal detection and the fractional backgrounds, as a function of and 
Qt, for Ufj, charged and neutral currents. Tables H, ^ and ^ summarize the results 
obtained (the cuts are the same than for the /i"*" analysis). Finally, Fig. ^ shows the 
signal efficiency and the fractional backgrounds for /i+'s and function of the 

neutrino energy. 

In summary, our study shows that a large magnetized iron calorimeter allows the 
detection, with high efficiency (~ 30-40%) of the golden-plated wrong-sign muon signal. 
The different backgrounds to this signal can be efficiently controlled using simple cuts, 
which exploit the different kinematics between signal and background events. The 
charm background can be suppressed to circa 10~^, taking advantage of the high degree 
of coUimation with the hadronic jet of charmed hadrons produced by 50 GeV muons. 
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Figure 10: Signal and backgrounds for CC and NC events. 



Baseline (km) 


CC 


77, CC 


77eNC 


77^ (signal) 


732 


6.9 X 10^ 


3.0 X 10^ 


3.1 X 10^ 


5 X 10^ 


3500 


3.0 X 10^ 


1.3 X 10^ 


1.4 X 10^ 


1.6 X 10^ 


7332 


6.9 X 10^ 


3.0 X 10^ 


3.1 X 10^ 


0.2 X 10^ 



Table 5: Data samples expected in a 40 kT detector for 10 useful jj, decays. 
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Ve CC 


+ Ve NC 


z/^ (signal) 


3.0 X 10-^ 


6.0 X 10-^ 


2.0 X 10~6 


4.0 X 10-1 



Table 6: Fractional backgrounds and signal selection efficiency for the wrong-sign muon 
search with fi" decays. 



Baseline (km) 


z/^ CC 


CC 


Uf, + Ve NC 


17^ (signal) 


732 


21 


18 


66 


2.0 X 10^ 


3500 


0.9 


0.8 


2.4 


6.4 X 10^ 


7332 


0.2 


0.2 


0.6 


8.0 X 10^ 



Table 7: Events surviving the cuts in a 40 kT detector for 10 useful /i decays. 

Instead, decays of energetic pions or kaons in NC events contaminate the signal at the 
10~^ level. The optimization of the S/N ratio points to the intermediate distances of 
0(3500) km as the optimal baseline. 



5 Analysis in energy bins 



A conservative estimate for the neutrino energy resolution in a detector of the type 
described in the previous section is AE^/E^, ~ 20%. For a /i beam of 50 GeV and the 
statistics of oscillated neutrinos expected in the range of parameters considered, it is 
reasonable to bin the data in five bins of equal width AE^ = 10 GeV. 

Let be the total number of wrong-sign muons detected when the factory is run 
in polarity p = /i"*",/!", grouped in 5 energy bins specified by i = 1 to 5, and three 
possible distances, A = 1 (732 km), 2 (3500 km), 3 (7332 km). 

In order to simulate a typical experimental situation we generate a set of "data" 
n^p as follows: for a given value of the oscillation parameters, the expected number of 
events, N^p, is computed; taking into account backgrounds and detection efficiencies 
per bin, b^ p and e^^, as given in Fig. |TT|, we then perform a Gaussian (or Poisson, 
depending on the number of events) smearing to mimic the statistical uncertainty. All 
in all, 

A Smear(Arr^e,^p + bf ) - b^p 
^i,P = ■ ^^^) 

H,p 

The "data" are then fitted to the theoretical expectation as a function of the neu- 
trino parameters under study, using a minimization. 
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Figure 11: Signal efficiency and total fractional backgrounds for the polarities /x"' 
and ii~ (down) as a function of the neutrino energy. 
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Figure 12: Asymptotic sensitivity to sin^ ^13 as a function of Am23 at 90% CL for 
L = 732 km (dashed lines), 3500 km (solid lines) and 7332 km (dotted lines), in the 
SMA-MSW solution. Only stastistical errors are included. 

where is the error of n^^ (we include no error in the efficiencies). We perform and 
compare six different fits using: x^, x\-, xi for the three distances, and the combinations 
Xi + xi' X2 + Xs! Xi + X2 + X3 to illustrate the gain in case the neutrino factory shoots 
to more than one location, a natural scenario given the ring configurations under study. 
For simplicity, we will consider a fit in at most two parameters at a time. All numerical 
results below will be obtained with the exact formulae for the oscillation probabilities. 



6 SMA-MSW or Vacuum solar deficit 



For the SMA-MSW or VO scenarios, the influence of solar parameters on the neutrino 
factory signals will be negligible^ and CP-violation out of reach. Besides its capability 



to ~ 1% m 



to reduce the errors on ^23 and 
machine to constrain/measure 6*13 

Consider first Oi-^. In Fig. |12|, we show the exclusion plot at 90% CL, on the Am; 



the factory would still be a unique 
and the sign of Ami, 12^ 



pg, 24 



(in the range allowed by SuperK) versus sin ^13 plane, obtained with the full unbinned 
statistics and the two polarities. The same results, but including as well background 



■^In practice, for the numerical results of this section, the central values in the SMA-MSW range 
are taken: Ato?2 = 6 x 10"^ eV^ and sin^ 26*12 = 0.006. 
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Figure 13: As in Fig. including as well background errors and detection efficiencies. 



errors and detection efficiencies are shown in Fig. The statistical treatment in the 
presence of backgrounds is done as in Notice that the sensitivity is better at L = 
3500 km than at 732 km when efficiencies and backgrounds are included. The latter 
are responsible for it. The sensitivity at 7332 km is also worse than at 3500 km, due to 
the loss in statistics. In conclusion, the sensitivity to di^ can be improved by three-four 
orders of magnitude with respect to the present limits. This is consistent with the 



results of [18| given the different statistics used. 



The second major topic would be to perform the ffist precise measurements related 
to matter effects, in order to determine the sign of Am23 ||2^ and the size of the matter 
parameter, A (xrie- 

We have studied the determination of the sign of Am23, assuming that the absolute 



value has by then been measured with a precision of 10%. 
region around the best fit values of SuperKamiokande: 



We have explored the 
= 2.8 X 10-3 eV^ and 



6'23 = 45°. We perform a x analysis on the Am23,6'i3 plane, as described in last 



section. The conclusion is that, for "data" generated within the range = 1-10° and 
|Am23| in the range allowed by SuperKamiokande, a missidentification of the sign of 
Am23 can be excluded at 99% CL at 3500 km and 7332 km, but not at the shortest 
distance, 732 km. This conclusion agrees with the analysis of ref. [^, which did not 
include the energy dependence information. 

We have further studied how the matter parameter of Eq. (O) and the angle ^13 can 



be measured simultaneously. Fig. 14 shows the result of a x fit as described in section 
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Figure 14: 68.5, 90, 99 % CL resulting from a simultaneous fit 0/613 ond A. The 
parameters used to generate the "data" are denoted by a star, while the haseline(s) 
used in the fit is indicated in each plot. Only statistical errors included. 



^. Only statistical errors have been included in this figure. The corresponding results 
including as well background errors and detection efficiencies are shown in Fig. 
At 732 km there is no sensitivity to the matter term, as expected. However, already 
at 3500 km, A can be measured with a 10% precision. At the largest baseline, the 
precision in A improves although at the expense of loosing precision in ^13 due to the 
loss in statistics. The level of precision discussed here might even be interesting for 
geophysicists 

The above conclusions also hold for the vacuum solution to the solar deficit. 



7 LMA-MSW 

Assume now the LMA-MSW scenario. Fixed values of the atmospheric parameters are 
used in this section, /S.m\^ = 2.8 x 10"'^ eV^ and maximal mixing, 623 = 45°. A precision 
of 1% in these parameters is achievable through muon disappearance measurements at 
the neutrino factory |^ . This level of uncertainty is not expected to affect the results 
of this section. 

Let us start discussing the measurement of the CP phase 6 versus 613. We have 
studied numerically how to disentangle them in the range 1-10° and < S < 180°. 

Consider first the upper solar mass range allowed by the LMA-MSW solution: 
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Figure 15: The same as in Fig. [I^ hut including backgrounds and efficiencies. 



Am^2 — 10~^ eV^. Fig. shows the confidence level contours for a simultaneous fit of 
6^13 and 5, for "data" corresponding to 6*13 = 8°, 5 = 54°, including only statistical errors 
in the analysis. Fig. [1^ shows the same analysis taking into account the background 
errors and detection efficiencies of Fig. 0. The correlation between 5 and ^13 is very 
large at the shortest baseline 732 km, as argued in section ^ The phase 5 is then 
not measurable and this indetermination induces a rather large error on the angle 
^13. However, at the intermediate baseline of 3500 km the two parameters can be 
disentangled and measured. At the largest baseline, the sensitivity to 5 is lost and the 
precision in di^ becomes worse due to the smaller statistics. The combination of the 
results for 3500 km with that for any one of the other distances improves the fit, but 
not in a dramatic way. Just one detector placed at O (3000 km) may be sufficient: a 
precision of few tenths of degree is attained for and of a few tens of degrees for 5. 

Similar figures are obtained for "data" corresponding to smaller values of ^13, as 
shown in Fig. |T8| for ^13 = 2°. The pattern is maintained as well for different values 
of 5: see Fig. 19 for "data" corresponding to 5 = 0° and ^13 = 8°. This last figure 
also proves that, if the sign of Am^2 is known by the time the neutrino factory will 
be operative, 5 = 0° is distinguishable from 6 = 180° with just one baseline. This 
exemplifies the power of the analysis of the energy dependence. Recall in any case that 
a TT-ambiguity on 6 has no bearing on the existence of CP-violation, and we will go on 
considering positive values of Am^2- 

The sensitivity to CP- violation decreases linearly with Am^2- the central value 
allowed by the LMA-MSW solution, Amj^ = 5 x 10"^ eV^, CP-violation can still be 
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7.6 7.8 8 8.2 7.6 7.8 8 8,2 



Figure 16: 68.5, 90, 99 % CL contours resulting from a fit of 9is and S. The 

parameters used to generate the "data" are depicted by a star and the baseline(s) which 
is used for the fit indicated in each plot. Only statistical errors are included. 
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Figure 18: As Fig. ( ^F^ for different values of the parameters denoted by the star. 
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Figure 19: As Fig. j^PJj ) for different values of the parameters denoted by the star. 
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Figure 20: Fit for Am^2 = 5 x 10 ^ eV^ including backgrounds errors and detection 
efficiencies. The star indicates the parameters used to generate the "data". 

discovered, as shown in Fig. At the lower value allowed, Am^2 = 1 x 10^^ eV^, the 
sensitivity to CP-violation is lost with the experimental set-up considered, as shown in 
Fig.0 

We have quantified what is the minimum value of Am^2 for which a maximal CP- 
odd phase, 6 = 90°, can be distinguished at 99% CL from S = 0°. The result is shown 
in Fig. |2^: Amf2 > 2 x 10~^ eV^, with very small dependence on ^13, in the range 
considered. 

One word of caution is pertinent: up to now we assumed |Am^2l sin 2^12 known 
by the time the neutrino factory will be operational. Otherwise, the correlation of these 
parameters with ^13 would be even more problematic than that between 6 and 613, as 
illustrated in Fig. ^ for |Am^2l- The error induced on the measurement of ^13 by the 
present uncertainty in |Amf2l is much larger than that stemming from the uncertainty 
on 6. Fortunately, LBL reactor experiments will measure |Amf2| and sin 2^12 if it lies 
in the LMA-MSW range. Even if the error in these measurements is as large as 50%, 
the problem would be much less serious. We have checked that such uncertainty does 
not affect our results concerning the sensitivity to 6, and only induces an error in 613 



that can be read from Fig. |23. 



Concerning the measurement of the matter parameter, we have considered simul- 
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Figure 21: As Fig. |^ for Amj^ = 1 x 10"^ eV^. 



L=3500 



a 9 10 



Figure 22: Lower limit in Am^2 which a maximal CP phase (90°) can be distinguished 
from a vanishing phase at 99% CL, as a function of Oi^ at L = 3500 km. Background 
errors and efficiencies are included. 



30 



Figure 23: Simultaneous fit to ^13 and Amfg- range shown in the vertical axis 

is the presently allowed LMA-MSW range. The star indicates the parameters used to 
generate the "data" and the CP-odd phase is set to zero. Backgrounds and detection 
efficiencies are included. 

taneous fits of ^13 and A, for two values of the CP-phase: 6 = 0,7r/2. The confidence 
level contours obtained are very similar to those in Fig. |1^ for the SMA-MSW solution. 
This indicates that there is no dangerous correlation between A and 613 in the presence 
of sizeable (5-dependent terms, and both parameters can be safely measured at 3500 
km. However, it is important to stress that the simultaneous measurement of the three 
parameters di^, 6 and A will increase the errors with respect to the two-parameter fits 
performed here. In this respect the combination of two baselines: (9(3500 km) and 
(9(7332 km) may be helpful. 
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8 Summary 



The neutrino beams obtained from muon storage rings will be excellent for precision 
neutrino physics. The appearance of wrong-sign muons is a powerful neutrino oscilla- 
tion signal, which allows to improve considerably our knowledge of the leptonic flavour 
sector. 

Two very important questions are the optimal energy for the decaying muons and 
the optimal detection distance(s), in view of the physics goals. The higher the parent 
muon energy, the larger the oscillation signals. This fact, together with the requirement 
of low backgrounds and good detection efficiencies, lead to consider muon energies as 
high as possible within realistic machine designs. Energies of a several tens of GeV are 
currently under discussion, assumed here to be = 50 GeV, for definiteness. 

Energy and detection distance are intertwined in the oscillation pattern of neutrinos 
propagating in matter. We have derived an analytical approximate formula for the 
oscillation probabilities in matter, which helps to understand how the sensitivity to 
the most interesting quantities scales with the neutrino energy and distance. 

We have shown that an analysis in neutrino energy bins, combined with a compar- 
ison of the signals obtained with the two polarities, allows to disentangle the unknown 
parameters at long enough baselines. In particular, for the LMA-MSW solution, 613 
and 6 can be simultaneously measured. We have also studied realistic backgrounds 
and detection efficiencies. The overall conclusion is that the intermediate baseline of 



0(3000 km) is optimal for the physics goals considered in this paper (see Fig. |2J for 
an artistic view of possible locations). 

Quantitatively, our two parameter fits at 3500 km indicate: 



• The angle ^13 can be measured with a precision of tenths of degrees, down to 
values of 6*13 = 1°. The asymptotic sensitivity to sin^6'i3 can be improved by 
three orders of magnitude or more. 

• If the solar deficit corresponds to solar parameters in the LMA-MSW range, 
CP-violation may be tackled. The phase 5 can be determined with a precision 
of tens of degrees, for the central values allowed for lAm^gl' maximal CP- 
violation can be disentangled from no CP-violation at 99% CL for values of 
|Am22| > 2 X 10-5 eV^. 

• The sign of the atmospheric mass difference, Am23, ^le determined at 99% 
CL, for ^^13 within the range ^13 = 1-10° and |Am23| range allowed by 
SuperKamiokande data. 

• A model independent confirmation of the MSW effect will be feasible, and the 
matter parameter A measured within a 10% precision, or better if combined with 
the longest baseline: 7332 km. 
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In the case of the LMA-MSW solution, the combination of the two longest basehnes 
may be useful if a multiparameter fit becomes necessary. 

Even though non-zero neutrino masses are barely established, the neutrino sector 
of the theory can be convincingly argued to herald physics well beyond the standard 
model. It is in this perspective that a neutrino factory should be built. 
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Figure 24: Artistic view of possible baselines of 0(3000 km) at which a neutrino factory, 
located at Geneva (Switzerland) or Geneva (Illinois), could shoot. The inner and outer 
circles correspond to 2000 and 4000 km, respectively. Note that the circles have not 
been projected onto the Earth surface. o . 



Appendix A: Non-oscillated statistics 

Tab. D contains the results for the Ve and fluxes for the parent muon energy 
= 50 GeV and for = 2 x 10^° useful muons per year, per 5 operational years. 
The result for three muon polarizations are shown: V^± = 0, =F 0.3 (the "natural" 



polarization) and =f1. For = our results agree with when the same number 
of useful muons is considered. 

For a /i+ beam of = 50 GeV, the average neutrino and antineutrino energies 
are, for V^+{^^~) = 0, {Ep^^^^)) = 35 GeV, and (^^..(p,)) = 30 GeV. For = 
-1(+1), we get = 30 GeV, = 30 GeV. 



= 50 GeV 






L (km) 


$.,/io- 




$p,/io- 


$../10^^ 


732 









123 


122 


123 


122 


±0.3 


110 


159 


110 


159 


±1 


81.5 


244 


81.5 


244 


3500 









5.35 


5.35 


5.35 


5.35 


±0.3 


4.83 


6.95 


4.83 


6.95 


±1 


3.56 


10.7 


3.56 


10.7 


7332 




P,^ 





1.22 


1.22 


1.22 


1.22 


±0.3 


1.10 


1.58 


1.10 


1.58 


±1 


0.81 


2.43 


0.81 


2.43 



Table 8: Neutrino and antineutrino fluxes for L = 732, 3500 and 1332 Km 'per m? per 
5 operational years when 2 x 10^° muons decay in the straight section of the storage 
ring. This fluxes have been averaged over an angular divergence of 0. 1 mr. 

Tab. 1^ contains the numerical results for the number of CC interaction rates for 
and /x^ fluxes in a beam in a 40 kT detector with = 2 x lO^'' useful muons 
per year per 5 operational years. These results represent the number of leptons of a 
given flavour observed at the detector in case no neutrino oscillation occurs, neglecting 
detection efficiencies. 
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= 50 GeV 




/^^ 


L (km) 










732 









692 


300 


352 


590 


±0.3 


603 


390 


306 


768 


±1 


394 


600 


200 


1180 


3500 









30.4 


13.1 


15.4 


25.8 


±0.3 


26.4 


17.1 


13.4 


33.6 


±1 


17.2 


26.2 


8.75 


51.6 


7332 









6.90 


3.00 


3.50 


5.88 


±0.3 


6.01 


3.88 


3.05 


7.66 


±1 


3.92 


5.98 


1.99 


11.8 



Table 9: Neutrino and antineutrino charged currents interaction rates for L = 732, 
3500 and 1332 km 'per 40 kT and per 5 operational years when 2 x 10^° muons decay 
in the straight section of the storage ring. These fluxes have been averaged over an 
angular divergence of 0. 1 mr. 



Appendix B: Oscillated statistics 

As an illustration, we give the oscillated fluxes for three values of the atmospheric 
mass difference, Am23 = 2, 4, 6 x 10~^ eV^. The rest of the leptonic parameters are 
Am?2 = 10"^ eV^, 623 = 45°, ^13 = 13° and Ou = 22.5°. The matter parameter is 
taken to be A = 1.1 x 10"^ eVVGeV for the baselines of L = 732 km and L = 3500 
km, while A = 1.5 x 10~^ eV^/GeV for L = 7332 km. The in principle measurable 
quantities are the number of leptons of a given flavour and charge reaching the detector 
at a given baseline. 



Tab. |l^ shows the total number of leptons of different flavours (z/g —>■ /i^ 



z/g — e , i/g — > — and i/^ — z/,- — > r for a /i beam) at 732 km. Tab. [TT] shows 
the analogous results from /i+ decays. Detection efficiencies are not included. The 



whole exercise is repeated for the baselines L = 3500 km and L = 7332 km in Tabs. [12 



13 and 14, RB, respectively. 
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Am^3 


iV^-/105 


iVe+/105 


iV^+/103 


iVe-/103 




iV^-/10^ 





0.002 


691 


300 


13.8 


24.1 


14.3 


20.6 


0.004 


684 


299 


51.3 


91.7 


52.9 


80.9 


0.006 


673 


298 


110 


201 


113 


178 


0.3 


0.002 


601 


390 


18.0 


22.0 


18.6 


19.4 


0.004 


595 


389 


66.7 


83.8 


68.8 


76.1 


0.006 


584 


387 


143 


184 


147 


167 


1 


0.002 


391 


600 


27.7 


18.6 


28.6 


16.5 


0.004 


387 


598 


103 


70.8 


106 


64.2 


0.006 


378 


596 


220 


155 


226 


140 



Table 10: Calculated charged currents event rates for jjL~ beam assuming neutrino os- 
cillations in a 40 kT detector, for a L — 732 km baseline as a function of E^j^, Amlg, 
for different polarizations of the parent muon. We have considered 1 x lO^-*" negative 
muons decays (2 x 10^° useful muons/year x 5 operational years). 





Amia 




Ne-/W' 






^"^-710^ 


A^r+/104 





0.002 


351 


590 


28.9 


11.7 


28.0 


10.5 


0.004 


348 


588 


110 


43.6 


106 


41.4 


0.006 


342 


586 


239 


94.1 


232 


91.5 


-0.3 


0.002 


305 


767 


37.5 


11.2 


36.4 


9.86 


0.004 


302 


765 


142 


41.9 


138 


38.7 


0.006 


297 


762 


311 


90.1 


302 


85.1 


-1 


0.002 


200 


1180 


57.8 


9.52 


56.0 


8.36 


0.004 


196 


1170 


219 


35.2 


212 


32.7 


0.006 


192 


1172 


478 


75.3 


464 


71.6 



Table 11: Calculated charged currents event rates for beam assuming neutrino os- 
cillations in a 40 kT detector, for a L — 732 km baseline as a function of E^, Am^g; 
for different polarizations of the parent muon. We have considered 1 x 10^^ positive 
muons decays (2 x 10^° useful muons/year x 5 operational years). 
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Am^3 




iVe+/104 


iV^+/102 


iVe-/103 




iV^-/103 





0.002 


282 


130 


65.8 


21.1 


75.3 


187 


0.004 


232 


128 


163 


81.9 


187 


631 


0.006 


170 


126 


227 


174 


259 


1150 


0.3 


0.002 


244 


169 


85.5 


17.9 


97.9 


176 


0.004 


198 


166 


212 


70.5 


243 


585 


0.006 


143 


164 


295 


150 


337 


1054 


1 


0.002 


156 


260 


132 


15.3 


151 


145 


0.004 


120 


256 


327 


58.7 


373 


463 


0.006 


81.0 


253 


454 


121 


518 


792 



Table 12: Calculated charged currents event rates for jjL~ beam assuming neutrino os- 
cillations in a 40 kT detector, for a L — 3500 km baseline as a function of Ef^, ^"^23; 
for different polarizations of the parent muon. We have considered 1 x lO^-*" negative 
muons decays (2 x 10^° useful muons/year x 5 operational years). 







iV^+/10^ 








iV^-/10''^ 


iV,+ /103 





0.002 


143 


254 


25.7 


60.4 


22.9 


100 


0.004 


118 


240 


96.4 


164 


88.1 


348 


0.006 


86.0 


221 


195 


249 


181 


655 


-0.3 


0.002 


124 


330 


33.4 


62.4 


29.8 


92.7 


0.004 


100 


312 


125 


165 


115 


319 


0.006 


72.4 


287 


254 


246 


235 


592 


-1 


0.002 


79.3 


507 


51.4 


50.2 


45.8 


77.0 


0.004 


60.8 


480 


193 


124 


176 


254 


0.006 


40.8 


441 


390 


173 


362 


450 



Table 13: Calculated charged currents event rates for jj,^ beam assuming neutrino os- 
cillations in a 40 kT detector, for a L = 3500 km baseline as a function of Ef^, Arri^^, 
for different polarizations of the parent muon. We have considered 1 x 10^^ positive 
muons decays (2 x 10^*^ useful muons/year x 5 operational years). 
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Am^3 


iV^-/103 


iVe+/103 


iV^+/102 


iVe-/103 









0.002 


527 


299 


2.38 


3.97 


1.94 


159 


0.004 


256 


296 


16.7 


28.2 


15.1 


406 


0.006 


89.1 


292 


35.8 


78.8 


33.5 


523 


0.3 


0.002 


451 


388 


3.10 


3.95 


2.52 


146 


0.004 


212 


385 


21.6 


27.3 


19.6 


362 


0.006 


75.1 


380 


46.6 


74.1 


43.6 


452 


1 


0.002 


273 


597 


4.77 


3.89 


3.88 


115 


0.004 


110 


592 


33.3 


25.2 


30.1 


257 


0.006 


42.5 


584 


71.6 


63.3 


67.0 


287 



Table 14: Calculated charged currents event rates for jjL~ beam assuming neutrino os- 
cillations in a 40 kT detector, for a L — 7332 km baseline as a function of Ef^, ^"^23; 
for different polarizations of the parent muon. We have considered 1 x lO^-*" negative 
muons decays (2 x 10^° useful muons/year x 5 operational years). 







iV^+/103 








iV^-/10''^ 


iV,+ /103 





0.002 


265 


577 


6.06 


1.16 


5.84 


85.6 


0.004 


125 


513 


38.1 


11.5 


37.8 


225 


0.006 


52.8 


399 


94.8 


31.2 


95.0 


295 


-0.3 


0.002 


227 


750 


7.88 


1.20 


7.59 


78.6 


0.004 


104 


667 


49.5 


11.0 


49.1 


201 


0.006 


47.6 


519 


123 


28.5 


123 


255 


-1 


0.002 


137 


1150 


12.1 


129 


11.7 


62.4 


0.004 


54.8 


1030 


76.2 


10.0 


75.6 


144 


0.006 


35.7 


798 


190 


22.3 


190 


161 



Table 15: Calculated charged currents event rates for jj,^ beam assuming neutrino os- 
cillations in a 40 kT detector, for a L = 7332 km baseline as a function of Ef^, Arri^^, 
for different polarizations of the parent muon. We have considered 1 x 10^^ positive 
muons decays (2 x 10^*^ useful muons/year x 5 operational years). 
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Appendix C: Perturbative Expansion of Oscilla- 
tion Probabilities 

In this Appendix, we describe the perturbative expansion we have performed to 
obtain Eq. (|T6|). The problem is to diagonahze the neutrino mass matrix, 



M = U 



/O 






Ai2 
























0/ 



(20) 






\0 Ai3, 

with U as given in Eq. (|1]). The exact diagonahzation of this matrix has been done in 
pO| . Since we are interested only in the case in which Am^2 is small, it is adequate 
to use perturbation theory to compute corrections to first order in this quantity. This 
leads to much simpler analytical formulae. 

0, the diagonahzation of this matrix is very simple P5 |, 




^MiA+B^ 



In the limit Am^2 



M- 



(0) 



/ Ai3±A-B^ 
' 2 



V 



f211 



The matrix of eigenvectors is 



where 9m^ is defined by: 



tan26'MT = 



Ai3 sin 26^13 



f22) 



(23) 



Ai3 cos 2^13 

9m^ is to be taken in the first (second) quadrant if A13 cos 26*13 =F^ is positive (negative). 

At first order in A12, the perturbation to Eq. (^Tj) (in the basis of non-perturbated 
eigenvectors) is: 

/O 0\ 

A12 U^U^. (24) 
\0 0/ 



The eigenvalues at first order in A 12 are: 

AS" 



(1) 



\(0) , 



^ Aiacos^^ 
+ C12A12 , 



M: 



3 - /V3°^ + S12 A12 sin^ 9m^ , 



Xl'> = A 

The corresponding (not normalized) eigenvectors are 



(25) 



,(1) 



^2 



,(1) 



sin 26*12 cos 6*^/^ A 



12 



'12 



Ai2sin 9M^e 



Ai3 + A - B, 



2B^ 



( sin 26*12 cos 9 A 



12 



sin 2^ 



12 



sin 6*^/^6 



12 



-Ai3-A + 5^ 

-sj2sin2^Ai2e±''5 



Ai3 + A + 5, 



2B 



- sin 26*12 sin ^M:p Ai2e 
A31 + A + fi. 



±i5 



(26) 



40 



where 0m^ = On - Om^- 



With these results, it is easy to compute the oscillation probabilities, keeping con- 
sistently terms up to first order in A12. We obtain: 



Py.u^{y.y^) = ^23 sin2(2^M^) sin^ ^ ^ 



2 2 

■^23 -^12 



sin(4^A. )sin(2^A. )sin2 



sin(26'i2) sin(26'23) sin(2^A/ ) sin 



n2 



sin2(2^j\r ) cos(2^a/:f) sin(5^ L) 



A12L" 



+ 



Ai2 X 



sm 



cos ±6 — 



^ sin sin 9m:^ 



(0) 



,(0) 



sin 6m^ sin 6m^ cos S sin ^ ^ 



A 



(0) 



(27) 



„2 „• 2/ofl \ ■„2 I ^ 



Pv^Ur(v,V.r) = C23 SlU (26'M:p ) SlU 



2 2 
C23"5l2 



sin(4^M^) sin(2^M^) sin^ ( ^) ^ 

B^L 



+ sin2(2^M^) cos(2^M^) sin(5^ L) 



A12L' 



sin(2^i2) sin(2^23) sin(2^Af:p) sin (^^— ) A12 x 

cos ±5 



sm 



A^^lA / Xf^ L\ f cos6]\i^cos6m^ sin 6 sin 6 



(0) 



' B L\ 1 
sin Oau sin cos 5 sin ' ^ 



A 



(0) 



,(0) 



(2^ 



1 - sin2(2^Afc)sin2 



B^L 



+ 



.2 

■'12 



sin(46'Af^) sin(26'Af:^) sin 



2,B^L^ A12 
5 



+ sin2(2^Af:p) cos(2^a/t) sin(fi^ L) 



(29) 



These formulae are valid for all values of ^13 and to first order in A12. It is rather 
straightforward to check that they reduce to the vacuum result for A ^ 0. 

In section 3, we have further considered an expansion in which not only A12 but 
also 6^13 are small. We have kept terms up to second order: i.e. (9(A^2^i3°)) C^(^i26'i3) 
and (9(A526'i3^). The latter two can be obtained from Eqs. ( P7| , pSD and (P5|), by per- 
forming an expansion in 6*13. On the other hand, the terms of (9(A^2^i3°) absent 
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in that approximation, they appear at next order in the expansion. They can be easily 
obtained, though, starting directly from the diagonalization of the mass matrix with 
^13 = 0: 

.±A+A^ \ 

M = U23{623)Uu{9'm^)\ ±A+A,,+C^ q ^/l2(^M^)^^/23(M^ (30) 

V Ais/ 



with = JAf2 + A'^T 2AAi2 cos26'i2 and 



• oa' - ^12 sin 2^12 

sm2^A'/^ = — . (31) 

The expansion of the corresponding probabilities to second order in A12 gives, 

2 2 n /^12\^ 2 fAL 

Pu.u^iu.u^) = C23sm 2^12 (^^j sm l^— 

2 2 n /^12\^ 2 fAL\ 

Pu^UriPeUr) = ^23 SIH 2^12 — j SlU — j 

P....(^ePe) = l-sin2 2^i2(^) sin^^^) . (32) 
From the first of these equations we obtain the term in of Eq. (p^G]). 
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